APS/123-QED 



Inverse freezing in the Hopfield Fermionic Ising Spin Glass 

S. G. MagalhaesE and C. V. Morais 
Departamento de Fisica, Universidade Federal de Santa Maria, 97105-900 Santa Maria, RS, Brazil 



on 
o 
o 

o 

Q 

On 
(N 

CO 



c3 



i 

c 

o 
o 



> 
ON 

ON 

o 



F. M. Zimmer 

Departamento de Fisica, Universidade do Estado de Santa Catarina, 89223-100, Joinville, SC, Brazil 

(Dated: December 29, 2009) 

In this work it is studied the Hopfield fermionic spin glass model which allows interpolating from 
trivial randomness to a highly frustrated regime. Therefore, it is possible to investigate whether 
or not frustration is an essential ingredient which would allow this magnetic disordered model 
to present naturally inverse freezing by comparing the two limits, trivial randomness and highly 
frustrated regime and how different levels of frustration could affect such unconventional phase 
transition. The problem is expressed in the path integral formalism where the spin operators are 
represented by bilinear combinations of Grassmann variables. The Grand Canonical Potential is 
obtained within the static approximation and one-step replica symmetry breaking scheme. As a 
result, phase diagrams temperature versus the chemical potential are obtained for several levels of 
frustration. Particularly, when the level of frustration is diminished, the reentrance related to the 
inverse freezing is gradually suppressed. 

PACS numbers: 



I. INTRODUCTION 

There is current interest (see, for instance, Refs. 
HllliiiS) in studying inverse transitions, melting or 
freezing, motivated by the quite unconventional situation 
present in such transitions in which the ordered phase 
is more entropic than the disordered one£. This appar- 
ent counter-intuitive transition has become even more 
interesting since there are now various physical systems 
displaying this kind of transition, as for example, mag- 
netic films 9 and, particularly, high-T c superconductors 10 . 
Therefore, the knowledge about which conditions are nec- 
essary for the existence of inverse transitions is a chal- 
lenging issueii. 

Among the proposed classical magnetic models which 
can show inverse transitions, the Blume Capel (BC)i and 
the Gathak-Sherrington (GS) 1 "^ models can be useful to 
clarify what would be such conditions. For instance, the 
phase diagram of the BC models displays inverse melt- 
ing with a phase transition from the ferromagnetic (FM) 
phase to the paramagnetic {PM) one when the temper- 
ature is decreased. However, in order to obtain inverse 
melting in the BC model, it is strictly necessary to impose 
in the problem the entropic advantage of the interacting 
state through parameter r — l/k > 1, where I and k are 
the degeneracy of ±1 and spins states, respectively^. In 
contrast with the BC model, the GS models - which is 
similar to the the BC model except that it has the spins 
couplings given as random Gaussian variables - shows 
naturally inverse freezing. In other words, there is a reen- 
trant first order boundary phase separating the spin glass 
(SG) and PM phases. However, there is no need of any 
entropic advantage of interacting states^. The previ- 
ous discussion suggests that the randomness is respon- 
sible for producing naturally inverse freezing in the GS 
model as compared with the inverse melting case in the 



BC model. In fact, to be more precise, it seems that the 
presence of non-trivial form of randomness, which means 
frustration, 14 would be the essential ingredient which al- 
lows inverse freezing in the GS model with no need of 
entropic advantage. That sets some questions: What is 
the actual role of non-trivial randomness to produce nat- 
urally inverse freezing? Is this natural inverse freezing 
robust when the level of frustration is diminished? Can 
models with trivial randomness present naturally inverse 
transitions? 

The purpose of the present work is to investigate what 
is the actual role of frustration as a basic condition to pro- 
duce naturally inverse freezing. That could be achieved 
by comparing two distinct situations, the trivial random- 
ness and the highly frustrated regime and also study- 
ing how natural inverse freezing would be affected when 
the level of frustration is varied. Therefore, the prob- 
lem is not only to find a random magnetic model which 
presents naturally inverse freezing, but also to get one 
which allows the level of frustration to vary. Quite re- 
cently, the fermionic Ising spin-glass (FISG) modeli has 
been proposed as a model able to display naturally in- 
verse freezing in a phase diagram temperature versus the 
chemical potential fi. This model has the spin opera- 
tors Sf given as bilinear combinations of creation and 
destruction fermionic operators which have four eigen- 
states, two of them non-magneti o 15 ' 16 : 17 . The spin-spin 
coupling Jij , likewise the GS model, is a random variable 
which follows a Gaussian distribution. Indeed, the exis- 
tence of inverse freezing in this model could be expected 
since there is a close relationship between the GS and the 
FISG models 18 . For instance, the partition function of 
these two models can be related by a mapping between 
the anisotropy constant D of the GS model and jj*^. 

The FISG model is also useful to examine how quan- 
tum effects, included by the presence of a transverse field 
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r 16 i 20 can affect inverse freezing. In Ref. 0, it has been 
shown that, when T is increased, the reentrance in the 
SG/PM first order boundary phase in the phase diagram 
temperature T versus the chemical potential /i gradu- 
ally disappears. This scenario also suggests that V in 
the FISG model plays the opposite role of the rather ar- 
tificial parameter r in the BC model discussed in Ref. 
1. Nevertheless, the FISG model, in its original form 
with Gaussian random couplings, is a strongly frustrated 
model. Therefore, it is not possible to investigate in this 
model how the change of the level of frustration could 
affect the inverse freezing. 

An alternative route to accomplish the previously men- 
tioned investigation would be to use the FISG model in a 
version in which one could adjust the level of frustration. 
In that sense, the classical Hopfield spin glas o 21 ' 22 can be 
quite useful. In this model, the interactions Jy between 
classical Ising spins are given below: 

where £f — ±1 (i or j — 1,2, ...,N, N is the number 
of sites) are independent random distributed variables. 
This type of interaction has been intensively used to 
study complex systems^. There are clearly two extreme 
cases for the classical Hopfield Ising spin glass model. 
When p = 1, it becomes the Mattis model 2 ^, which is a 
well known example of trivial randomness 1 —. The second 
one is when p — > oo. In that limit, the thermodynam- 
ics corresponds to the strongly frustrated regime given 
by a random Gaussian distributed Jif^- In particular, 
the replica symmetric mean field solution of the classical 
Hopfield spin glass model 2 ^ can provide a useful method 
for the purposes of the present work. It is described in 
terms of two order parameters, the usual spin glass q 
and Trig = 53i=i £f < Si > which indicates the pres- 
ence of Mattis states that have a thermodynamics similar 
to the usual ferromagnetis m 21 i 22 . Most important, there 
is a parameter called here degree of frustration, a = p/N 
which allows controlling the level of frustration in the 
theory. It is important to remark that there is a particu- 
lar value a c which for a > a c the effects of frustration are 
dominan t 21 ' 22 . As a consequence, it is possible to inter- 
polate the thermodynamics from the trivial randomness 
(a = 0) to the strongly frustrated regime (a > a c ). 

Therefore, to answer the questions arisen previously, 
we use the Hopfield Fermionic Ising Spin Glass (HFISG) 
model. In this case, the spin operators Sf are defined 
as the FISG model, but the random spin coupling is 
given as Eq. {!]). The partition function is obtained in 
the functional integral formalism using Grassmann fields 
where the disorder is treated with the replica method 25 . 
The problem can be reduced to a one-site problem by 
using a similar procedure to solve the classical Hopfield 
spin glas a 21 ' 22 within the static approximation (SA) 26 
and one-step replica symmetry breaking (1S-RSB)2!. It 
should be remarked that in the HFISG model, as the 



FISG on o 15 i 17 ' 20 i 28 , the replica diagonal component of 
the SG order parameter q aa appears as an additional or- 
der parameter to be solved with the others. Particularly, 
this order parameter represents the spin self-interaction 
that has an imaginary time dependence for quantum SG 
model o 26 ' 29 . However, the spin operators Sf commute 
with the Hamiltonian operators for the present HFISG 
model. Consequently, the q aa has no dynamic, which 
means that the SA is exact in the present wor k 15 ' 16 . 

One important point in the present work is how to lo- 
cate first order boundaries phases found for different de- 
grees of frustration. The criterion adop ted here follows 
closely that one suggested in Ref. yO for the classical 
GS model. Hence, we selected, from the set of spin glass 
solutions in the transition, that one which meets contin- 
uously with the spin glass solution for small value of the 
chemical potential /i. In fact, that corresponds to the 
largest spin glass order parameter which gives the low- 
est Grand Canonical Potential (see also the discussion in 
Ref. l3lh . Then, by equating the Grand Canonical Po- 
tential of the spin glass and paramagnetic solutions, the 
first order boundary phase is located. 

The use of the 1S-RSB scheme also deserves some re- 
marks since our main interest is to obtain boundaries 
of phase transitions. It has been found in previous 
works 2,5,7 that the use of the replica breaking symmetry 
(RSB) schemes in GS or FISG models essentially pre- 
serves the reentrance in the SG/PM first order bound- 
ary phase associated to the inverse freezing. There are 
only small differences which appear mainly at very low 
temperatures. Even so, for the HFISG model, we de- 
cided to use 1S-RSB scheme to check for intcrmedianted 
values of a whether or not these differences will remain 
unimportant. Moreover, we also analyzed the stability of 
the RS solutions for a sake of completeness of the work. 

This paper is structured as follow: in Section 2, we de- 
rived the thermodynamics and the set of coupled equa- 
tions for the saddle point order parameters. In Section 
3, phase diagrams temperature versus the chemical po- 
tential are presented for several values of a. The entropy 
behavior as function of temperature and the grand canon- 
ical potential as function of chemical potential are also 
discussed. Finally, Section 4 is reserved to conclusions. 



II. MODEL 

The Hamiltonian considered here is a Hopfield FISG 
(HFISG) model 

H = -Y,JijS!Sl (2) 

ij 

where Jij is given in Eq. {T]) and the random £f follows 
the distribution 

m) = \ + \ (3) 
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In this model, Sf = h[fii\ — HjjJ is the spin operator, 
with fiia — cla-Cia as the number operator, cJ CT (ci a ) are 
fcrmions creation (destruction) operators and a =f or j 
indicate the spin projections. 

The partition function in the grand canonical ensemble 
is given in the Lagrangian path integral formalism where 
the spin operators are represented as bilinear combina- 
tions of anticommuting Grassmann fields (</>, (j>*) 32 



Z{fi} = J D( 



where 



(4) 



(5) 



H(P(t),<Kt))} 



where /i is the chemical potential and (3 = 1/T. Z{fi} is 
Fourier transformed in time which results in 



Z{ M } = / D{4>*<t>)e 



with 



(6) 



(7) 



A sg = J2H PJiiSnJ^i-w'), (8) 

Q ij 



and 



where /i is chemical potential, ui — (2m + 1)tt, lo = 
2mir (m = 0, ±1, ■••) are the Matsubara's frequencies 
and a s = +(— ) if =T (I) • In this work, the problem is 
analyzed within SA which considers only the term when 
J = in Eq. 

The Grand Canonical Potential is obtained by using 
the replica method: 

(3n = -]iml/(nN)(((Z(n))) ( -l) (10) 

n— >0 

where Z(n) = Z n and means the configurational 

averaged over £. Thus: 

Z(n) = J] / exp[^ + A%%*\ (11) 

0=1 ^ 

where after using given in Eq. (Q]), the action Agg* 
can be written as 



si c 



SG 



/3J 

2iV 



EE(E^*) 5 



£=1 a— 1 



2iV 



(12) 



with a denoting the replica index and Sf = (0). 

The average over Z(n) given in Eq. (JTTJ) is discussed 
in detail in the Appendix A. In the present work, the 
one-step replica symmetry breaking (1S-RSB) ansatz is 
adopted, in which the replica matrix {Q} and the matrix 
{r} are parametrized as: 



q if a = (3 
c/„ . \ qi if |a - /3| < x, r aP 
qo otherwise 



r if a = j3 
r\ if |a — j3\ < x 
ro otherwise 

(13) 

and order parameters m" are invariant with respect to 
permutations of replicas: mf — m, where a — 1, • • -n. 
Therefore, the parametrization (fT3")) is used in Eqs. (|A8[) - 
(|A9p . As a consequence, the 1S-RSB Grand Canonical 
Potential is found as 



(30, = -(3(i - 



(3 2 J 2 a 



(rq - (1 - x)rxqi - xr q ) 



a ln 1 ~ PJ[q - gi + x(gi - g )] 
2x n \-pJ(q- qi ) 



1 



f3Jaq 



(3Jm 2 



21- f3J[q- qi +x{qi-qo)] 

1 



a 

+ 2 



(14) 



x ln[l - 0J(q - ft)] - lim - ln((9({r}, m, 

n— »0 n 



where 



s ?(<"') = JeE^L^+^m. (9) x exp 



e({r},m,0 = 
f3 2 J 2 a 



l//3J)E(S a ) 2 



x exp 



x exp 



n/x 

(ri -r )E 



E 



a=l 



5*" 



1=1 \a=(l-l)x+l 
n n 

r (E^) 2 +^E(^ Q ) 



(15) 



with a = f>/AT. The quadratic forms into the function 
0({r},m,£) can be linearized by Hubbard-Stratonovich 
transformations where new auxiliary fields are introduced 
in the problem. Therefore, one has: 



e({r},m,f) = / Dz 



Dv{ Dw 



D{4?4>) expEC^)^ 1 ^)^^) 



x-i n/x (16) 



with Dy = dy ^ 27r (y — z, v, w) and 

G-^u) = g-\u) + h(z, v, w) + (3J(£m). (17) 
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The local spin glass component of the random field 
h{z,v,w) is defined by 



h(z, v, id) — 0J[y / aroz + y/ a{r\ — tq)v 



(18) 



(r — n - 1/0 J)w]. 

The functional integral over the Grassmann variables, 
as well as the sum over the Matsubara's frequencies in Eq. 
(I16[) . can be performed following closely the procedure 
given in reference o 16 ' 20 . Finally, the Grand Canonical 
Potential within 1S-RSB approximation can be written 
as 

0Jm 2 



i 



(1 - x)nqi - xr q ) + 
0Jaq Q 



2 1 - 
-Qi)\ 



0J[q 
a 



2x 



1 



In 



a 

- qi + x(qi - g )] 2 
1 - 0J[q-qi + x(qi - 



ln[l - 0J(q 
-Qo)] 



(19) 



l-0J(q- qi ) 



0V~- I Dz((\n / Dv(2K{z,v\0T))i 



where 
K{z, v\£) = cosh/3/i - 

and 



e Ma^J(r-rx)-l] 



coshH (z, v, £) 
(20) 



H(z, v\$ = 0J[y/E^z + y/a(n - r )v + £m]. (21) 

The set of equations for the order parameters m, qo, qi, 
q and the block size parameter x can be found from Eq. 
(fT9"]) using the saddle point conditions. 

In particular, the elements of matrix r are given by 

qo 



{1 - 0J[q - qi + x( qi - qoW 



(22) 



n - r 



qi - qo 



- 0J{q - qi )]{l - 0J[q - qi + x( qi - q )}} 

(23) 

r ~ ri = 0J[l-0J(q-qi)Y 

The average over £ in the Grand Canonical Potential 
can be done using the parity properties of the functions 
dependent on z and v. The entropy can also be obtained 
directly from the Grand Canonical potential. 

The stability of the RS solution q a p = q, r a p = r 
(a 7^ 0) is studied using the Almeida-Thouless analysis, 
in which the condition for the stable RS solution ( X a t > 
0) is obtained as^ 



X at - [1 + 0J(q - q)f 

e u cosh H(r) 



a(0J)' 



Dz 



K(r, r) 



e u sinhff(r) 
K(r, f) 



(25) 



with K(r, r) — cosh(/3/i) + e" cosh H(r) and 



u = -0Ja[0J(f-r)-l], H(r) = 0Jy/E?z+0Jm. (26) 



III. RESULTS 

The numerical solutions for the coupled set for the sad- 
dle point order parameters q%, qo, q and m are displayed 
in phase diagrams T/J versus \x/ J given below for sev- 
eral values of the parameter a, where T is the temper- 
ature and J is defined in Eq. (fT]). For the numerical 
results, J = 1 is used. The RS scheme for the order 
parameters is also calculated when qo = qi = q. The sta- 
bility of such solution is investigated by calculating the 
Almeida-Thouless eigenvalue Xat given in Eqs. (|25|) - 
(f2"6")l . Such analysis shows that the PM solutions and 
Mattis states (FE) present Xat > while the SG ones 
have Xat < 0. This result could indicate that it is neces- 
sary to use RSB schemes to locate more adequately the 
SG/PM first order boundary phase. However, the re- 
sults shown in Figs. (U)-® have found that, in such first 
order boundary phase, the reentrance and, therefore, the 
inverse freezing are not essentially affected by the use of 
RS or 1S-RSB schemes similar to the FISG modeli. 

In Fig. (TT|) , it is presented phase diagrams which illus- 
trate two quite distinct situations concerning the degree 
of frustration. The first one (a = 0), which is shown in 
Fig. ((T]-a), corresponds to trivial randomness 14 , which 
means that there are no effects of frustration. In this 
phase diagram, for low T/J and small fx/J, one gets 
FE as solution. For higher T/J and/or larger fi/J, it 
is found the PM solution. For small /x/J, the boundary 
phase transition between FE and PM phase (called here 
T 2c (/x)) is second order. However, when /x/J increases, 
it appears a tricritical point located at T tc = J/3 and 
Htc = 0.438 J (see Appendix B). Most important, there 
is no reentrance in the subsequent first order part of the 
boundary phase transition Ti c ([i). In Fig. (JTJ-b), the de- 
gree of frustration is increased (a = 0.1) and, for that 
case, effects of frustration become dominant. Thus, the 
FE solution is replaced by a SG one in which, m = 0, 
q 7^ 0, qo 7^ and q\ ^ 0, while the PM phase has 
m = 0, q 0, go = and qi — 0. The freezing tem- 
perature has also a second order part T 2 /(^). As for the 
trivial randomness tricritical point appears, with 

value T tc ~ 0.438 J and fitc ~ 0.752 J (see Appendix B). 
Nevertheless, the first order part of the freezing temper- 
ature T\ f (/i) now shows a reentrance which indicates the 
existence of inverse freezing, as it can be seen in the in- 
sert in Fig. (QJb), which displays the entropy as function 
of T / J. There, it is shown that entropy of the SG phase 
is larger than the PM one for that particular degree of 
frustration. It should be remarked that this phase di- 
agram is quite similar to that one found for the FISG 
model-. The location of the boundaries phases in Fig. 
(fT]) has been checked in the limits [i = and T = as it 
can be seen in Appendices B and C. 

There is a more complex scenario as compared with 
that one described in Fig. (QJb) when a — » 0. In Fig. 
J2]-a) (a — 0.01), the solution for the order parameters 
shows a phase diagram which illustrates this new sce- 
nario. For a small /x/J, when T/J is decreased, there is 
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FIG. 1: Phase Diagrams T/J versus /x/J for a = 0.0 and a = 0.1. In the phase diagrams, the full lines represent the second 
order transition and the dashed lines represent the first order transition. Ttc indicates the tricritical point. The insets show 
the entropy S versus T/J for specific values of fi/J in the first order transition, where T lc and T\f represent the first order 
transition between FE/PM and SG/PM phases, respectively. 
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FIG. 2: Phase diagrams T/J versus fi/J for several low values of a. T tc is a tricritical point and Tt p is a triple point. Dashed 
lines represent the first order transition and solid lines represent the second order transition. The insets show the location of 
the triple point in details. Close to the triple point, a small difference between the first order boundary of the 1S-RSB (dotted 
lines) and RS (dashed lines) solutions can be seen in the insets. 



a second order phase transition between the PM and SG 
phases. However, for even lower temperatures, there is 
another phase transition, which is now a first order one, 
between the SG phase and FE region, which is given now 
by m 7^ 0, g 7^ 0, 9o 7^ and q\ ^ 0. For this particular 
value of a, SG and FE solutions occupy approximately 
equal sizes in the phase diagram. The freezing temper- 



ature has a similar behavior to that one found in Fig. 
((T]-b). It has a second order part T 2 f(n) for small fJL, 
then it appears a tricritical point at Tt c ~ 0.366J and 
fitc ~ 0.539 J (see Appendix B). Below this point, Tif(fi) 
presents a reentrance which allows, for an adequate con- 
stant fi, crossing from the SG phase to the PM one when 
the temperature is decreasing. Nonetheless, this first or- 
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der boundary phase has a complex nature. It appears a 
triple point at (T tp ,fi tp ) where FE, PM and SG phases 
coexist. Below this point, the first order boundary phase 
Ti c (/i) diplays no reentrance, as in the case a = 0. 

For smaller values of a, as shown in Figs. ([2]-b)-([2]- 
d), the region where FE solutions are found becomes 
increasingly larger than the SG one which is consistent 
with earlier results found in the classical Hopfield spin 
glass model which displays, in a phase diagram T ver- 
sus a, a dominance of the Mattis states when a — » O^i. 
Even so, in this new situation, the location T^f{p) is 
not affected so much. However, the tricritical point is 
displaced for smaller and lower values of fi/J and T/J, 
respectively. In comparison, the triple point is displaced 
for smaller values of \x/J and higher values of T/J. As 
a consequence, Tif(/j,) and T lc (/x) appear in a decreasing 
and increasing range of temperature, respectively. Nev- 
ertheless, most important, the reentrance in Xi/((/i) is 
gradually suppressed when a — ► 0. However, the insert 
in Fig (121-d) shows that even when the level of frustration 
is very weak, and T2f(fi) appears in a very short range 
of temperature, a reentrance in such first order boundary 
phase is still preserved. 
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FIG. 3: Grand canonical potential Q versus fi/ J for T/J — 0.2 
and two values of a = 0.1. The full, dotted and dashed- 
dotted lines represent the £1 of SG, PM and FE solutions, 
respectively. 
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The labels Ti c (p) and Ti/(/x) indicate the first order transition 
between the FE/PM and SG/PM phases, respectively. The 
full vertical line indicates the second order transition between 
SG/PM phases. 



which, as in the inset of Fig. ((T]-b), indicates the exis- 
tence of inverse freezing. The same procedure has been 
adopted in Figs. |[23-b)-(|23-d), values of \i are adjusted to 
be close to \i tv . The result found in Fig. (|2]-b) is similar 
to that one found in Fig.([2]-a). However, Fig. (|2]-c) dis- 
plays the entropy behavior when T\f{p) and T lc (/z) are 
now crossed. The first crossing is in the reentrant line 
transition T\f(fi) giving an inverse freezing as in Figs. 
d2J-a)-([2]-b). The second one in Ti c (/z) gives a usual phase 
transition in which the PM phase is more entropic than 
the FE one, as already found in Fig. (Q~|-a). ln Fig- (0- 
d) the reentrance in Tif(p) is almost suppressed. Thus, 
for the chosen value of fi, there is only one crossing in 
Ti c (/x) which gives the entropy behavior of a usual phase 
transition as that one found in Fig. JT]-a). 



The procedure to locate the first order boundary lines 
in the previous phase diagrams is illustrated in fig. (3), 
where the grand canonical potential versus \i/ J is plotted 
for T/J = 0.2. In fig. (3), the set of multiple SG solu- 
tions is represented by branches of full lines, where the 
chosen SG solution is that one which meets continuously 
with the only one SG solution available for small values 
of jit/J. From that SG solution, the first order bound- 
ary phase is obtained by equating the grand canonical 
potential of SG, PM and FE solutions. 

The corresponding behavior of the entropy as a func- 
tion of T/J for the values of a used in Fig. ([2]) is shown in 
Fig. (j4|). This figure illustrates the gradual suppression 
of inverse freezing when a — * 0. For instance, in Fig. ^ 
a), the value of /i = 0.53 J is chosen to cross the reentrant 
first order boundary phase in Ti/(/x) in Fig. ((2]-a). The 
figure shows that the entropy of SG phase is larger than 
the entropy of the PM one in the first order transition 



IV. CONCLUSION 

In the present work, it has been studied the HFISG 
model in which the random spin-spin coupling Jjj, in- 
stead of the usual Gaussian distribution, is given by Eq. 
(Q]) . For this particular choice of Jy , the problem can be 
treated in a mean field framework which allows adjust- 
ing the level of frustration. Therefore, it is possible to 
study not only the existence of a natural inverse freezing 
in the limits of trivial randomness and strong frustration 
but also how changing the level of frustration could af- 
fect such transition. In this approach, within 1S-RSB 
scheme, besides the replica non-diagonal SG orders pa- 
rameters go j Qi aud the block size parameter x, there is 
also the replica diagonal order parameter q. The set of 
order parameters is completed with m which corresponds 
to the presence of Mattis states (FE) 21 i 22 . The coupled 
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equations for q , qi, x, q and to are solved for several 
values of degree of frustration a in a phase diagram T 
versus chemical potential fi. In particular, the RS solu- 
tion is also obtained when qo = q\ = q. 

The comparison among the several scenarios displayed 
in Figs. CE}-© elucidates the role of frustration as the 
essential ingredient responsible for producing naturally 
inverse freezing in the HFISG model. In other words, an 
inverse freezing without any need of an artificial entropic 
advantage. For instance, the comparison between results 
given in Figs. {T]-a) and (JT]-b) shows that only in the case 
a = 0.1 (the strong frustrated regime) a reentrant first 
order boundary phase appears. The trivial randomness 
is not able to generate naturally a reentrance and, there- 
fore, an inverse transition. Furthermore, in Fig. I[2|). 
when the level of frustration is diminished, the range in 
temperature of the SG/PM first order boundary phase 
is decreased but, most importantly, the reentrance in this 
first order boundary phase is gradually suppressed, what 
implies in the gradual suppression of inverse freezing un- 
til the complete disappearing in the trivial randomness 
case (a — 0). Nevertheless, whatever the level of frus- 
tration, there is always a reentrance in SG/PM first or- 
der boundary phase which gives an inverse freezing. In 
contrast, the range of the FE/PM first order boundary 
phase is simply increased, nothing else happens. These 
previous features lead to the conclusion that frustration 
in any level is the necessary condition to create naturally 
the entropic advantage of the SG phase as compared with 
the PM one which generates an inverse freezing. Albeit, 
these results are restrict to a particular model, we sug- 
gest that the role of frustration as an essential ingredient 
to rise naturally inverse freezing could be more general. 

One last remark must be done. The present approach 
could be used directly in GS model with the coupling 
Jij as given in Eq. ([1]) replacing the original Gaussian 
distributed onei£. However, the great advantage of the 
HFISG model is that it would also allow studying the 
present problem in the presence of a transverse magnetic 
field Di&aa^. In Ref. 0, it has been shown that T tends 
to destroy the inverse freezing. Therefore, the presence 
of an additional T would lead to another question: what 
would happen with the inverse freezing in the HFISG 
model when the level of frustration and the strength of 
quantum effects are simultaneously changed? In that 
case, it would possible to investigate the robustness of 
the inverse freezing by adjusting simultaneously the de- 
gree of frustration and T. This question is currently un- 
der investigation. 
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APPENDIX A: THE AVERAGE OVER Z(n) 

In this appendix, the averaging procedure of the parti- 
tion function of the HFISG model is introduced following 
closely Ref. |2j. The first term in the action A s qq (see 
Eq. (TT2|) ) can be linearized by a Hubbard-Stratonovich 
transformation by introducing n x p auxiliary fields to" 
which are splitted in two subsets with nx (p — l) and nxl 
terms. Therefore, 



exp(A^ 4 ) = cxp 



2N 



EEm 2 



/CO ( j 

-°° I u=l a 

poo f P 

J —OQ ii-i 



-,«) 2 + vl 



g=l+l a 



1 



where r 



(Al) 



din',] 



hY,iC e) Sfm<* v{e) and 



Dm "(e) ~ n„( e ) IL ^/fr- 
it is assumed that the relevant contributions come from 
to" which are order unity, while to^ is of order 1/y/N. 
Therefore, the average over the p — l independent random 
variables £f can be done using P(£f) given in Eq. ([3]) 
which results in: 



((exp 



g=l+l a i 



exp 



EE ln (<x*b.(PJ y Z,S?m°) 



(A2) 



The argument of the exponential in the right hand side 
of Eq. (|A2|) can be expanded up to second order in . 
The result is a quadratic term of the spins variables Sf 
in the last exponential of Eq. (|A1[) . This term can be 
linearized by introducing the spin glass order parameter 
q a /3 using the integral representation of the delta function 
as 
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dr 



a/3 



2tt 



exp 



(A3) 



Therefore, the exponential involving to" in Eq. ()A1[) 
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can be written as: 
exp | pN £ E 

I g=l+l a 



5K) 2 + ^(E^K 



n 

a/3 



dq a pdf af: 



2vr 



exp 



9 £ m^A^m^D^ 



e=i+i 



with 



a/3 i 

where the matrix element 

A Q /3 = (1 - f3q aa )6 al3 + 0q a p(l - 5 a p) . (A6) 

Introducing Eqs. (|A"4"1) - (|A"6)) into Eq. (|aT) . the to^ fields 
can be integrated to give: 



((exp(^ t ))) c =exp 



2N 



/+oo ( 1 
-oo , i „, 



"2«) +<« 



n 

a/3 



^|^ex P {^-^-0TrlnA)}. 



(A7) 



Assuming / = 1 in Eq. (|A7[) . the averaged partition 
function is given as 



iiry, NU /°° n a Z" 30 TT dq a pdf a p yj dq aa df aa 
((Z(n)))z= I Dm, I [[ p - [[ - 



a^/3 



x exp 



x exp 



' ^ ^ ^aaQaot ~t~ ^ ^ ^ ^"a^Qa/: 



/3JA 



a#/3 



p-1 



TrlnA 



((0(^,r ao ,m?)»^ 



(A8) 



where the matrix A is defined in Eq. (|A6|1 and 
9(r a p ) f aa ,mZ)=f[ ( D(<t>* a ct> a )exp[A? 

+ E E E ^« Ms^M'^w (<<• 



(A9) 



with 

A? 



a#/3 



E( 5 f ) 2 + ^ EE x 



(A10) 



and c/- = + /3fi. The trace of the matrix A given 



in Eq. (1A6|) is obtained in terms of its eigenvalues. 

The Grand Canonical Potential is found introducing 
Eq. (|A8[) m Eq. (fT0|) which is evaluated at the saddle 
point. Thus, 



—5— <(«*?) ) = — ^~ 



(All) 



(A4) an d 



r a /3 



— — ((mfmf)) = 



pr Q/3 ; a^. (A12) 



APPENDIX B: TRICRITICAL POINTS 

This appendix presents a procedure to obtain the tri- 
critical point (^t c ,T tc ) for both cases: a — and a > 0. 
For a = 0, the Landau expansion of the Grand Canonical 
Potential, Eq. (|19[) . in powers of m gives 



j3Vt(m) = A Q + A 2 m 2 + A^m A + A 6 m 6 + 



(Bl) 



with 



A 2 = (3J[1 - (/3J)/(1 + cosh(/3/z))]/2, (B2) 

A 4 = /3 4 J 4 [2 - cosh(/3^)]sec/i[(/3/i)/2] 4 /96. (B3) 

The tricritical point is obtained when A 2 = A4 = 0. The 
position of T 2c in the Fig. (Q~|-a) for fi = can be also 
checked from the equation A 2 = 0. 

For a > 0, the PM/SG phase transition is investi- 
gated. In this case, it is assumed that there is no essential 
difference between RS and 1S-RSB schemes concerning 
the location of the tricritical point. 2 Therefore, we start 
with the thermodynamic potential within the RS solution 
(q = qi = qo, r = ri = ro and x = 0) written explicitly 
as a function of the order parameters r and f, and with 
m = 0: 



f3Q = -pn + %\8Jf In (f - r)l - Dz 

2 r — r J 

x ln[cosh/?jU + e'^ -5 ( r - r )-^2~\ coshf3Jy/arz] 



(B4) 



where r and f are given by Eqs. (|22| and (|24|) within the 
RS solution, respectively. 

Equation (|B4[) is expressed as an expansion in powers 
of r, which is related to the SG order parameter. There- 
fore, 



(B5) 



where f(r, /z, T, a) is obtained by a saddle point solution 
of f2. In this case, r(r, n, T, a) can also be written in the 
form of a series 



r = Tq + T\T + T2T 



(B6) 
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with n = as a result, 



and 



Xn 



f = [^(i-wr 1 , 

2 + /3 6 J 6 a 2 f % 2 (l~X ) 
>o[2-/3 4 J 4 af 2 X (l-X )]' 



exp 



") 



e /3Ja/2 CQsh ^ + exp ( /3 2 J-W 



(B7) 
(B8) 

(B9) 



Now, Eq. (|B6j) is introduced into the coefficients of 
expression (|B5I) , which are expanded in powers of r again. 
The resulting expression is then expressed in powers of 
the SG order parameter q by expanding r: 



pn = f q + 



with 



F 2 q 2 



F 3 q 3 



F 4 q 4 



(1 - l3Jq ) 4 (1 - PJqo) 6 (1 - PJqo) 8 

(BIO) 



F 2 = > 4 J 4 a * 2 *) 



/;,,: = | (/ ?WX 3 + ^) 



(BH) 
(B12) 



For a = 0, the grand canonical potential of FM solution 
at T — is fi/ m = —fi — 4f. By comparing fi/ m with the 
grand potential potential of PM phase, fl pm = —2/i, the 
first order boundary is located at fj, — 0.5 J. For the cases 
shown in Fig. it is recovered the grand canonical 

potential for a = 0. 

For a > a c , the effects of frustration are dominants. In 
this case for RS solution, q — q oc T for T — » 0. Therefore, 
close to the transition, which means fj,/J > a P J (9-q) 
the fi sg at T = is 



2[l-/9J(g-g)] ■ 



y \ e^ 2 /2 

v er/c 1 TiJ " W 



(CI) 



where erfc(z) = 1 
2/ = 7^ (m/^- 



-erf(z) (erf(z) is the error function), 
and r = erfc (~^J ■ 



From the condition Q sg — fi pm , y and r can be solved 



0). Particularly, fJ,i/(T 



which allows finding fiif(T 
0)=0.689 for a = 0.1. 

In 1S-RSB scheme, for T -> 0, x = - ?i) is 

independent of T and x = jT/J. Therefore, close to the 
transition, where fi/J > a P J ^~i 1 ) 
is 



2[l-/3J(g- 9 i)] ■ 



the n sg at T = 



/?WX 2 



12 



(45X 2 - 12X + 1)} (B13) 
-i[(ror 2 - I) 2 + h - 6/3J(l - /8Jgb)F 3 

where g = + \/a)]. The PM/SG second order 

phase transition occurs when i*2 = and -F4 > with 
the tricritical point located when F 2 — and F4 = 0. In 
particular, one can use the condition F 2 — with Eqs. 
(B7) and (B9) to obtain the critical temperature T 2 f by 
solving 

cosh/3 2/ M = exp ( J/32 2 /V ° )[J/Ml + V5) - 1] (B14) 

where (3 2 f = 1/T 2 f. For fi = 0, it is recovered the result 
T 2/ = 0.729 for a = 0.1 



APPENDIX C: FIRST ORDER TRANSITION AT 
ZERO TEMPERATURE 

Here a procedure is presented to obtain the first or- 
der boundary phases of FE/PM phases (a = 0) and 
SG/PM phases (a > a c ) at T = within the RS and 
1S-RSB solutions. 



fi.. 



Ja 

T 



9i 



1-X 



+ 7(7-1 gi - r qo) 



1-X - 7(9i - 9o) 7 



1 , /" 1 - X - 7(9i - 9o ) 



(i-x) 



where 




(C2) 



(C3) 



^2a(ri - r ) 



with ?y± = i^/ar^z + 07 (n - ro) — /i, /i = A 4 / J - f t=^: 
and $± = ±^/aroz+ §7(^1 — ro) — M- The previous results 
for RS solution are recovered from eqs. (|C2|) - (|C3p when 
q = qi = qo. Again, equations for x, 9o, <Zi, »"o, ri and 
7 can be obtained from Q sg and solved. In this case, 
fii f (T = 0)=0.682 for a = 0.1. 
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